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SOME EXPERIMENTS ON BATEMAN-HORN 
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Abstract. We describe some studies related to the frequency of 
prime values of integer polynomials. 


1. Introduction 


The Bateman-Horn conjecture (see m) states that given m irre¬ 
ducible polynomials with integer coefficients, then the num¬ 

ber of positive integers k < x such that fi{k) is prim^^ for every 
1 < z < m is asymptotic to 

c r dt 

Dj, 

where D is the product of the degrees of the fi, while 


c(/)=n 

p 


1 -^p(/)/p 

i-i/p 


where np{f) is the number of solutions of f{x) = 0 mod p, where 

f{x) = 

In this note, we look at the simplest case where m = 1, and study 
the behavior of the coefficient C{f), where / is a random monic poly¬ 
nomial in Z[x]. Our model of a random polynomial is one where all 
the coefficients (except the leading one, which is always equal to 1) 
are uniformly chosen from the interval [—A, A], where A is typically 
equal to 1000. Further, we approximate the constant C{f) by taking 
the products over the hrst several primes (hrst 3000 primes in the ex¬ 
periment below). Now, there is also a hnite probability that C{f) = 0 - 
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^We allow our primes to be negative. 

^Such a polynomial is almost surely irreducible, though in the experiments we 
throw away the few exceptional cases. 
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this would be true if the values of / are always divisible by p, for some 
prime p - obviously the most likely such prime is 2, so it makes sense 
to look at those / which don’t always vanish mod p for any p (or, at 
least, any p we look at)j^ 

Now, what do we look at? 

• What is the mean value of C{f) over our sample space? (whether 
all irreducible / or those / satisfying the Bunyakovsky condi¬ 
tion? 

• How are the values of C{f) distributed (here, it only makes 
sense to look at those / satisfying the Bunyakovsky condition). 

• What are the extremal values of C, and what property of the 
polynomials involved is responsible? (we have not studied this 
yet). 

It should be noted that in the “orthogonal” direction (where we have 
a number of linear polynomials, which corresponds to the prime fc-tuple 
conjecture), some statistical results were obtained by P. X. Gallagher 
[2] and E. Kowalski [3]. 

2. Some experimental observations 


I looked at polynomials of degrees 2 through 6. Here are the means: 


degree 

mean 

Bunyakovsky mean 

Bunyakovsky log mean 

2 

1.00329 

1.33525 

0.0847562 

3 

1.00714 

1.38343 

0.129013 

4 

0.971798 

1.35929 

0.124901 

5 

1.00044 

1.38399 

0.132433 

6 

1.02139 

1.40398 

0.155772 

7 

1.02448 

1.41654 

0.166858 

8 

1.00884 

1.36817 

0.120447 


It does not seem surprising that the mean of all the irreducible poly¬ 
nomials is very close to 1. However, I am unaware of any proof of this. 
Here is a somewhat relevant fact: 


Lemma 2.1. For each f G Fp[a;] let n{f) be the number of zeros of f 
counted without multiplicity. Then, the mean value of n{f) over all of 
Fp[x] is 1. 

Proof. Consider the variety V defined over ¥p by 

n—1 

X'^ + Y^ YiX^ = 0 . 

i=0 

^This is the so-called Bunyakovsky condition. 
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The statement of the Lemma is equivalent to asserting that there are 
p" Fp-points on V. However, if we rewrite the equation of V as 

n—1 

X^ + '^V,X^ = -Vo, 

i=l 

it becomes obvious that V is nothing but the affine space A"’, so the 
result follows. □ 

Now, in our experiment we sample uniformly from polynomials in 
'L[x\ with coefficients bounded above by C and below by —C, then 
reduce modulo p. When C is much bigger than P our sample space 
will contain all of Fp, but with somewhat uneven multiplicity. When 
C is much smaller than p, our sample space contains only a small 
sliver of Fp, so we are asserting that the Fp points of V are very well 
equidistributed. If we average over primes, such a result is true. 

Lemma 2.2. Let f G Z[a:] be irreducible over Q. Then, 

lim ^ 

x^-oo Tl[x) 

where np{f) is the number of zeros of f modulo p. 

Proof. For each prime, let Sf = {di,... ,dk} be the set of degrees of 
irreducible factors of / modulo p. Chebotarev’s theorem (see, e.g., m) 
says that the fraction of the primes for which a given s/ occurs is the 
same (asymptotically) as the fraction of the elements of the Galois 
group of / which have the cycle decomposition with lengths given by 
Sf. If f is irreducible, the Galois group of / is a transitive subgroup of 
Sn- The result now follows from Lemma [2. 3 1 □ 

Lemma 2.3. The average number of the fixed points of elements of a 
transitive permutation group eguals 1. 

Proof. This follows from Burnside’s lemma: for a group G acting on a 
set X, we have 

IMGI = 2;; Z A"l. 

where |W/G| is the number of orbits of the G action, and \X^\ is the 
number of fixed points of g. Since G acts transitively, the left hand side 
equals 1. □ 





4 


IGOR RIVIN 



(a) Histogram of log C{f) 



(b) quantile-quantile plot of log 

Cif) 


Figure 1. Degree 2 



(a) Histogram of log C{f) 



(b) quantile-quantile plot of log 

C{f) 


Figure 2. Degree 3 


2.1. The distribution. As you may or may not be convinced by the 
charts below, the values of C{f) for Bunyakovsky / (that is, an / 
which does not vanish identically modulo any prime p) seem to be log- 
normally distributee^ If true, this seems to indicate that the terms 
1 — Up/p are independent; the distributions for different p are not iden¬ 
tical, but are the ones given by Chebotarev density. Now, independence 
for (very) small primes follows from the Chinese Remainder Theorem, 
but when the primes are large compared to the coefficients of our poly¬ 
nomials, that is far from clear. 

3. How PRIME-RICH CAN A POLYNOMIAL BE? 

As we have found (at least experimentally), the average value of 
C{f) is 1, which means that a garden-variety monic (irreducible) poly¬ 
nomial of degree d, when given an argument of size around N, is about 
as likely to give a prime value as a number of size around N'^. We say 
that a polynomial is prime-rich if it is much more likely to give prime 

^The graphs are of the Bateman-Horn quantity C/D, not of C. 
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(a) Histogram of log C{f) 

Figure 3. 



(b) quantile-quantile plot of log 

C{f) 


Degree 4 



(a) Histogram of log C{f) 



(b) quantile-quantile plot of log 

C{f) 


Figure 4. Degree 5 



(a) Histogram of log C{f) 



(b) quantile-quantile plot of log 

Cif) 


Figure 5. Degree 6 


values (which means that C{f) >1). Of course, since just over a quar¬ 
ter of all polynomials fail the Bunyakovsky condition (the probability 
that a polynomial with content 1 in Z[a:] gives even values always is 
around 1/4, odd values always is 1/27, and so on), the average Bun¬ 
yakovsky polynomial is somewhat prime rich, but how well can we do? 
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(a) Histogram of log C{f) 

Figure 6. 



(b) quantile-quantile plot of log 

C{f) 


Degree 7 



(a) Histogram of log C{f) 

Figure 7. 



(b) quantile-quantile plot of log 

Cif) 


Degree 8 


To find the answer, we computed 5000 random polynomials of coeffi¬ 
cients bounded above by 5000 of each of the degrees 2, 3,4, 5, and found 
the ones most prime-rich. Below, we give the three richest for each of 
the degrees we computed, but not before asking: 

Question 3.1. Is C{f) bounded above (hxing degree)? 

Of course, if the distribution of C{f) is genuinely log-normal, then 
the answer is NO. 

One place to look for an answer is in the proof that the product 
dehning C{f) actually converges. The proof (see [Tj) uses Landau’s 
Prime Ideal Theorem to observe that 

= loglogx + Af + o{l), 

p<x ^ 
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for any irreducible / (the constant Af for /(x) = x is the so-called 
Mertens-Meissel constant M). In turn, that implies that 


hm y ^ A,- M, 

x^oo ^ p 

p<x ^ 


which immediately implies that the product converges, and indicates 
that the size of C{f) is controlled by the magic constant Aj, which 
does not seem to help that much. 

Note that if we had some uniform estimate on the convergence speed 
of the Euler product, then the answer to the Question 3N is clearly 
negative: simply take a polynomial which has no roots for the hrst k 
primes (this is easy to do by Chinese remaindering) - the product of 
the hrst k terms of the dehning product is going to be of order of log k 
(while the coefficients of the polynomial with given reductions mod p 
will be of order of e^) - a uniform estimate on the remainder would 
hnish the job. However, we are not aware of any such uniformity result 
(notice that it would be enough to have uniformity on average). If 
we remove the condition of being monic, then the obvious candidate 
(without Chinese remaindering) is the polynomial 


Pn,k{.x) = JJpi. 

i=l 


A numerical experiment indicates that (^(Paog) ~ 9.5. Modulo all of 
the above, one can ask a more detailed question: 

Question 3.2. If we look at all / G Z[x] of some hxed degree with coef- 
hcients bounded (in absolute value) by N, is it true that the maximal 
C{f) is of order log log iV? 

One can look further, and look at whether more sophisticated meth¬ 
ods of computing C{f) than just multiplying out the beginning of the 
Euler products give us any clue. Such methods are described in a 
number of very nice papers by Nobushige Kurokaw8[^ - [H CJ El [6] 
- Kurokawa shows how to represent C{f) as a product of Artin L- 
functions. Again, it is not clear how to leverage this to get a bound. 

The other obvious question is: 


Question 3.3. What do prime-rich polynomials have in common? 

It is fairly clear {assuming the truth of the Bateman-Horn conjecture 
) that the constant term should be prime (or a product of large primes), 
and this is borne out by the results below. 

Here are the winners (we give top three for each degree): 
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C{f) 

fix) 

6.3722 

- 2619a; + 1291 

6.36569 

a;2 - 2717a; - 1471 

6.23592 

a;2 + 2321a; + 911 

5.51225 

a;3 + 3914a;2 - 3485a; + 2773 

5.37671 

a;3 - 611a;2 - 424a; - 2999 

5.35378 

+ 707a;2 - 800a; - 509 

6.19895 

x^ - 1065a;3 + 409a;2 - 265a; + 817 

5.74642 

X* + 254a;3 - 3125a;2 - 1204a; - 2609 

5.65711 

x^ + 4590a;3 - 4932a;2 + 2061a; - 1289 

6.27693 

a;® - 2127a;^ - 1190a;3 - 2317a;2 - 1499a; - 2257 

6.15153 

x^ - 4686a;^ + 2812a;3 - 4475a;2 - 1714a; + 3389 

5.74885 

x^ - 3738x^ - 150a;3 - 4819a;2 - 2748a; + 307 

6.86439 

x^ + 2697x^ + 3377x^ + 2484a;3 - 2700a;2 + 1587a; + 1831 

6.28748 

a;® + 4705a;® + 232a;^ - 3661a;3 + 3063a;2 - 820a; - 3533 

6.12438 

a;® + 1138a;® - 2757a;^ - 4376a;® + 954a;2 + 4060a; - 2729 

5.72267 

a;^ + 1965a;® + 2378a;® - 2384a;^ - 1298a;® - 600a;2 - 2610a; - 1249 

5.62225 

a;^ - 4618a;® + 4170a;® + 4299a;^ + 1447a;® + 4695a;2 + 2032a; - 4387 

5.54845 

a;^ + 704a;® - 2286a;® - 1938a;^ - 462a;® + 2470a;2 + 4241a; - 2179 
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